Transient optical response of ultrafast nonequilibrium excited metals: Effects of 
electron-electron contribution to collisional absorption 
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Approaching energy coupling in laser irradiated metals, we point out the role of electron-electron 
collision as an efficient control factor for ultrafast optical absorption. The high degree of laser in- 
duced electron-ion nonequilibrium drives a complex absorption pattern with consequences on the 
transient optical properties. Consequently, high electronic temperatures determines largely the col- 
lision frequency and establish a transition between absorptive regimes in solid and plasma phases. 
In particular, taking into account umklapp electron-electron collisions, we performed hydrodynamic 
simulation of laser-matter interaction to calculate laser energy deposition during the electron-ion 
nonequilibrium stage and subsequent matter transformation phases. We observe strong correlations 
between optical and thermodynamic properties according to experimental situations. A suitable con- 
nection between solid and plasma regimes is chosen in accordance to models that describe behavior 
in extreme, asymptotic regimes. The proposed approach describes as well situations encountered in 
pump-probe type of experiments, where the state of matter is probed after initial excitation. Com- 
parison with experimental measurements shows simulation results which are sufficiently accurate 
to interpret the observed material behavior. A numerical probe is proposed to analyse transient 
optical properties of matter exposed to ultrashort pulsed laser irradiation at moderate and high in- 
tensities. Various thermodynamic states are assigned to the observed optical variation. Qualitative 
indications on the amount of energy coupled in the irradiated targets are obtained. 

PACS numbers: 52.38.Dx,42.25.Bs,05.70.Ln,79.20.Ds 



I. INTRODUCTION 



Laser-induced phase transitions are extensively investi- 
gated by means of optical measurements due to their fun- 
damental interest and potential in material sciences and 
engineering. Optical probing was identified as a relatively 
straightforward way to time-resolve changes of material 
properties. However, singular effects in the material re- 
sponse in extreme conditions are not always easy to iden- 
tify in the sequence of matter alterations [l], [§, d, 0, @ . 
Although more accurate measurements using X-ray or 
electron diffraction have emerged to probe matter trans- 
formation with atomic resolution 0, 0] , optical probing 
is still a powerful technique to resolve evolution changes 
in material transformations in a wide range of ther- 
modynamical states. First investigations of ultrashort- 
pulsed laser energy relaxation in metals concentrated on 
a wide range of effects including pump-probe reflection 
and transmission efficiencies Q, time-resolved particle 
emission and two-photon photoemission experiments Q . 
Subpicosecond lasers have been used to resolve nonequi- 
librium heating in Au and Cu [lol . [Tl| and to measure the 
electron-phonon coupling for a variety of metals [l]| [l3[ . 
Modelling the correct behavior of the material proper- 
ties under strong excitation is still a challenge for under- 
standing transport properties and relaxation processes in 
nonequilibrium systems [Til HH, [lfij] . Irradiated by an in- 
tense laser beam, a material undergoes a complex series 
of transformations, from solid-solid, solid-liquid, liquid- 
gas, up to the plasma phase. Each transformation has 



a particular optical signature and the description of the 
evolution remains a fundamental issue [13, [3, Ell l20| - 
Hence, comparison of experiments with dedicated mod- 
els is needed to validate the modelling approach even 
though experiments on extreme conditions of matter are 
difficult to perform and the data are scarce and often 
contradictory. Nevertheless, since experiments involve 
complex interconnected phenomena, various models with 
different degrees of complexity are required to decode the 
experimental information. 

A specificity of ultrafast laser excitation is its highly 
transient nature. Ultrafast laser irradiation can push 
matter in extreme high pressure, high temperature, 
nonequilibrium states where hot electrons interact with 
colder ions [2l|. Hereby, nonequilibrium describes the 
temperature mismatch between hot electrons and ions, 
while each interacting subsystem is considered thermal- 
ized on the timescales involved. If this is a crude approx- 
imation in the first moments of interaction, the estab- 
lishment of a high collision rate will reduce the domain 
of incertitude. Subsequently, the properties of the ma- 
terial are governed by the relaxation paths towards the 
equilibrium state [22|]. Moreover, since optical probing 
only involves the electronic system, information on lat- 
tice structural changes must be inferred indirectly from 
changes in the electronic system. From a fundamental 
and theoretical point of view, interesting features emerge 
when the pulse width becomes shorter than the electron- 
phonon relaxation time scale for energy deposition into 
the solid (23| . For this reason, we will be approaching sys- 
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terns whose dynamical evolution exhibits a pattern driven 
by the electron-electron and electron-ion collisional inter- 
actions (such as reflectivity or absorption coefficients). 
The energy distribution among each species determines 
changes within the excited material which affect the rel- 
ative transport properties. Since optical properties de- 
pend strongly on the ionization degree which is in turn 
controlled by temperature and density evolution, a hy- 
drodynamic approach is required. A coupled electromag- 
netic and hydrodynamic formalism, associated to specific 
optical and thermal models reflecting noncquilibrium fea- 
tures, provides an effective way to reproduce laser-matter 
experiments. Such simulations have been performed us- 
ing the ESTHER code, developed by the Commissariat 
a l'Energie Atomique, France. Details to the code as- 
sumptions and approaches are given in Ref. 24.25], and 
will be briefly reviewed in the next sections, focused on 
absorption and thermodynamic description. 

A full theoretical description of solid-laser interaction 
requires the selfconsistent treatment of hydrodynamics 
and the evolution of the transport properties which are 
associated. Despite many efforts, a global view including 
the multitude of interaction and relaxation processes is 
still incomplete and only particular aspects are under- 
stood at this time. The analysis of the transient optical 
properties of a laser-induced dense nonequilibrium plas- 
mas requires a description of the electron degenerecency 
effects. Different approaches have been used to model 
interaction mechanisms and subsequent energy absorp- 
tion. For considering the role of these interactions on 
degenerate matter properties, ionic and electronic struc- 
ture effects in relation to transport properties have to be 
calculated [III H3, [H, Hi| . Some of the studies dealing 
with ultrafast absorption phenomena are focused on the 
complex effect of simultaneous contribution of collision 
processes on the kinetic equation [l6|, [3(| HH . In consid- 
ering the nonequilibrium dynamics of electrons in met- 
als, the mutual influence of (e-e) and (e-p) interactions 
remains an issue of interest. Although quantum mechan- 
ical approaches and kinetic treatments have been devel- 
oped to caracterize the absorption rate in dense plas- 
mas and solids [32l l33l l34j and significant progress was 
made, the formulation of the electrons excitation through 
all these regimes remains difficult and limited. Conse- 
quently, we concentrate our calculations on the effects 
resulting from bringing together basic models to discrim- 
inate dominant mechanisms in several experimental sit- 
uations, encompassing a large variety of thermodynamic 
states. For excitation below certain relaxation times 
(electron-phonon equilibration), nonequilibrium appears, 
particularly in the femtosecond excitation regime. The 
dominant effects on the nonlinear optical dynamics ar- 
rive from the Pauli exclusion principle and the Coulomb 
electron-electron interactions among the free carriers ( e- 
e). The phonons also play an important role, especially 
on picosecond time scales. During and following their 
excitation, the electrons undergo a number of different 
relaxation stages as they scatter among themselves via 



the Coulomb interaction and with the phonons via the 
electron-phonon ( e-p ) interaction. In optical absorption, 
contribution may appear from electronic processes that 
allows momentum variation and conductivity changes. 
The ( e-e) interaction contributes indirectly to the inverse 
bremsstrahlung process of absorption since it screens the 
electron-ion pseudo-potential, and directly via umklapp 
processes. For simplicity and according to some other in- 
vestigations dedicated to optical probing [3, [H, [36| , we 
adopt a phenomenological approach as being well suited 
for integration in hydrodynamic codes. In this natural 
way, collision rates are described under the relaxation- 
time approximation for both (e-e) and (e-p) (37j . Obvi- 
ously, this simplified approach can not describe specific 
features of the collision process, but, along side with dif- 
ferent kind of experiments, it is adequate to improve ad 
hoc transport models in solid material up to the plasma 
phase. This will help reproducing transient evolution 
of the optical properties. By measuring the amplitude 
of the optical field emerging from the irradiated sample 
upon reflection, one can gain valuable insight into the 
interplay between many-body interactions and quantum 
confinement effects during very short time scales. 

The aim of this paper is to provide a deeper insight 
into the transport characteristics and their consequence 
on the optical properties in nonequilibrium systems, such 
as ultrashort laser-excited metals. With support of nu- 
merical calculations, we concentrate our investigations on 
several issues. First of all, the effect of electron-electron 
collision on the transient optical properties is modelled. 
Our original motivation was to insert a more accurate de- 
scription of the laser absorption in the ESTHER code. In 
order to validate our modelling assumptions, we perform 
simulations of reflectivity dependence on the angle of in- 
cidence, polarization and intensity while comparing the 
results with experimental data available in the literature 
[H, [39[ . Our results reveal that calculations match the 
experimental data with great accuracy when consider- 
ing the effect of (e-e) collision rate enhancement during 
laser irradiation. To extend our investigations, a great 
challenge is to compare the time evolution of the calcu- 
lated optical properties with time-resolved experiments. 
In fact, with subpicosecond laser pulses, it is also pos- 
sible to probe in a time-resolved fashion the charge car- 
rier dynamics in metals by a pump-probe optical system. 
Recent measurements on reflectivity and transmittivity 
evolution under ultrafast laser irradiation in pump-probe 
type of experiments have been reported ||. To realize 
a similar numerical pump- probe experiment, further de- 
velopments are required to identify potentially-induced 
phenomena in the framework of an optical, thermal and 
hydrodynamical study. Straightforward numerical mod- 
eling of such problems is a difficult challenge due to the 
various nature of underlying physical phenomena. While 
performing an analysis of the measured optical proper- 
ties, we propose an interpretation of the results by split- 
ting the behavior into two stages : a dense and nonequi- 
librium plasma phase following by a hot expanded plasma 
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regime. Our laser-matter numerical tool demonstrates 
new abilities to investigate nonequilibrium and associ- 
ated collisional effects on the material states ranging from 
cold solid to hot plasma. 

The rest of this paper is organized as follows. We firstly 
detail the modelling of the laser absorption by the mate- 
rial in Sees. II and III. In Sec. IV, we make use of com- 
parison with available experimental data sets to validate 
our numerical modelling. We then broaden the range 
of simulation by investigating the electromagnetic inter- 
actions relative to pump-probe numerical experiments. 
These experiments ranging from high intensities to be- 
low ablation threshold excitation, furnish arguments that 
demonstrate the applicability of such a numerical tool. 
Conclusions are drawn in the last section. 



II. ELECTROMAGNETIC EXCITATION BY 
ULTRASHORT LASER PULSE 

A. Laser-matter interaction background 

When electromagnetic irradiation is shorter than sev- 
eral ps, the equilibration time between the solid lattice 
and the laser-heated electrons is longer than the laser 
pulse length. Thus, thermal equilibrium could not be en- 
sured and the thermal state of the solid lattice is not well 
known. In order to understand the complexity of mate- 
rial absorption under the laser irradiation, the modelling 
of laser-matter interactions requires the interconnection 
of different physics fields, including optics, heat transfer, 
fluid dynamics, and phase transitions in nonequilibrium 
frames. Among different models that treat the interac- 
tion between ultrafast laser light and metallic targets at 
moderate intensities and provide scenarios of material 
ejection [l6|, HO, lU 5H , the hydrodynamic approach 
employing tabulated Bushman-Lomonosov-Fortov [iij 
multiphase equations of state (EOS) addresses reason- 
ably well the complexity of the problem. It also has the 
potential to reveal plausible evolution paths for the ex- 
cited matter that have also relevance for improving the 
laser interaction process [25L [45|. The ID hydrodynamic 
code ESTHER includes a Two- Temperatures Model to 
reproduce electronic thermal diffusion which takes place 
during electron-ion nonequilibrium. It is extended by a 
nonequilibrium hydrodynamic model decoupling electron 
and ion pressure components which both contribute to 
the ultrafast matter deformation according to the Euler 
equations. Fluid equations are solved to account for the 
conservation of energy, momentum, and mass in a La- 
grangian formalism. The specificity of this approach is 
the flexible dimension of the cells which deform together 
with the material and permit an accurate description of 
shock propagation. The simulation accounts for mechan- 
ical processes and particularly the effects on material 
strength that occur in shock propagation. This is due to 
constitutive models which describe the pressure and tem- 
perature dependence of shear modulus and yield strength 



[401 . This code is enhanced with a wave-equation solver 
to calculate the electromagnetic field in the material dur- 
ing laser-matter interaction. To determine optical re- 
sponse, especially in the presence of a density gradient, 
Maxwell equations are reduced to the Hclmholtz equation 
by using slowly-varying-envelope approximation for the 
electric field. This is justified because, for a laser pulse 
of r = 150 fs, the wavelength, A = 800 nm, is smaller 
than the coherence length L r = cr, where c is the speed 
of light . 



B. Optical property modelling 

If the idealized density step-profile is no longer actual, 
the Fresnel equations cannot be used and it becomes nec- 
essary to calculate the absorption in a non-uniform den- 
sity profile. The electric field is calculated numerically 
by solving the wave equation, using a one-dimensional 
mesh. The cell coordinates evolve in time due to the 
Lagrangian approach used to solve the hydrodynamic 
equations. For two pulses excitation, the contribution 
of each pulse is accounted for. The subsript a refers to 
the pump laser field at the frequency to a and (3 refers to 
the probe one at the frequency ujq. We consider a planar 
wave propagating along the z-axis. We denote z? the an- 
gle which the normal to the wave makes with the z-axis. 
Each component (cc, y, z) of the electric field is defined by 
the equation S a p — A^ y ^ z (z) exp(ifc° ^xsm.9) where 

a = u a,a/c is the vacuum wave number. We write the 
following Hclmholtz equation which is solved two times 
for the complex amplitudes of each components (a, 0) of 
the electric fields: 

Here, the complex vector k depends on the dielectric 
function e: 

~kl,3 = (k° a , f[e-sm^}. (2) 

In the solid case, the complex vector k is defined as a 
function of the refractive index rj and the extinction co- 
efficient k : 

k~a,/3 = ( k a,f3) 2 [(v^a ~ iK a ^) 2 - sin 2 i? Qlj a]. (3) 

r\ and k are taken at room temperature as a function of 
the wavelength from Ref. [47| . 

When the density decreases at the metal surface, the 
optical properties are strongly different from the above 
values and the electrical conductivity a for dense metal 
plasma at the equilibrium temperature is determined us- 
ing the tabulated values given in Ref. ;48] . Experimental 
measurements of the electrical conductivity of dense cop- 
per and aluminum plasmas in the equilibrium tempera- 
ture range of several eV and for a density range from the 
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solid density down to 0.01 g/cm 3 [4^| have been com- 
pared with theoretical models for conductivity by Ebel- 
ing et aLfUl]. We used the resultant ouc and the ion- 
ization rate under a tabulated form for given electron 
density and temperatures in our calculation. A phe- 
nomenological collision frequency is deduced from these 
values to determine the collision frequency at equilibrium 
v e q{N e , T e = Ti). It has to be noted that an interpolation 
is performed between solid and plasma phase for inter- 
mediate regime. For the liquid, gas and plasma cases, 
the dielectric function is given by the Drude model and 
the equation @ becomes: 



Here, Sq is the vacuum dielectric permittivity. 



(4) 



C. Reflectivity and transmissivity calculations 

The reflectivity 1Z P ' S , where (p,s) denotes the polari- 
sation state, is determined in a general manner, by cal- 
culating the ratio between the reflected and the incident 
power. The transmissivity ^? p > s is evaluated by the ratio 
between the transmitted and the incident power. The 
absorptivity is determined by the ratio between the ab- 
sorbed and the incident power. The relation 3~ v,s — 
1 — 1Z P,S — A p,s is satisfied. The angle dependence of 
the incident laser absorption is an important factor for a 
good optimization of the interaction. In this numerical 
experiment, the reflection factor is determined as the av- 
eraged value integrated on the measurement time. The 
calculation is realised by weightcning the incident inten- 
sity with the following form : 



IO*)l * 



(5) 



^» = J -^~r (6) 

As mentioned, simulation results are obtained in 
performing pump-probe numerical experiments during 
which the Helmholtz equations is solved for both pump 
and probe pulses. In this context, the reflectivity and 
the transmissivity can be associated to distinctive wave- 
lengths corresponding to each kind of pulse. In fact, ma- 
terial optical indices evolve as a function of the corre- 
sponding wavelength in time for each pulse. In the Es- 
ther code, both reflectivity and transmissivity are cal- 
culated at each time step St of the simulation. This al- 
lows that, according to the specific type of experiment, 



these parameters are either integrated over the duration 
of the pulse or they give momentary values. In the for- 
mer case, t = tl is used for the integral up boundary. 
For experiments providing data resolved in time below 
the picosecond timescale, r = St has been used in the 
previous expressions. 



D. Absorption collisional aspects in stages 
subsequent to excitation 

Electrons interact via the bare Coulomb interaction 
and after time intervals of the order of the inverse plas- 
mon frequency corresponding to the relevant carrier den- 
sity, the interactions become screened. Immediately after 
the photon absorption phase, collisions among electrons 
cause loss of coherence and a hot population of electrons 
is formed. Electrons and ions are initially described by 
two independent populations. The electrons evolve into 
hot thermal distributions characterized by an electron 
temperature that can exceed significantly that of the lat- 
tice. Such a time evolution, as the carriers equilibrate 
among themselves, is mainly induced by the (e-e) and 
(e-ph) interaction. The duration of this stage is deter- 
mined by the energy relaxation time, which can be of the 
order of picoseconds. The carrier dynamics affects the 
absorption properties through the conductivity parame- 
ter. In a first approximation, this quantity is completely 
described by the frequency-dependent Drude contribu- 
tion : 



N e e v + iu} a> p 
ml v 2 +u 2 a ^ 



(7) 



where the collision frequency v is given by the sum of 
the electron-electron (v ee ) and electron-ion [y el ) contri- 
butions. We can therefore express the collision frequency 
as a sum of equilibrium and noncquilibrium contributions 



y neq 



= v eq (N ei Ti) + i/ neq (N e ,T e ) - v neq (N e ,T e = Ti\8) 

Here, v eq (N e ,Ti) and v neg (N e ,T e = Ti) are (e-i) colli- 
sion frequencies at equilibrium, calculated in two differ- 
ent ways. The first term is determined using the tabu- 
lated values of conductivity at equilibrium [481 ] . These 
tabulated values were shown to match well experimental 
situations 49]. The last term corresponds to the equi- 
librium contribution to the collision frequency calculated 
by our nonequilibrium model. To ensure a convergence 
towards tabulated equilibrium values, we substract the 
equilibrium contribution from the nonequilibrium one. 
In this manner, the contribution of electrons out of equi- 
librium is treated as a correction of the known values of 
the equilibrium conductivity. Finally, the second term, 
l/ neij(A r e, T e ), is the nonequilibrium collision frequency for 
electrons which account for (e-e) (via umklapp processes) 
and (e-i) interactions for an electronic temperature larger 
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than the ionic one. This later contribution is calculated 
independently from the former one as it will be seen in 
the next paragraph. 

The electron-electron part is a fundamental issue for 
the description of low-temperature or nonequilibrium 
transport phenomena. It is a difficult problem to deal 
correctly with electron-electron scattering in calculations 
of transport coefficients. One of the possible way is the 
Boltzmann transport equation. However, the Boltzmann 
transport equation cannot usually be solved exactly and 
some approximations have to be introduced, such as a 
relaxation-time approximation. No tractable model ex- 
ists to describe the conductivity on a large range of tem- 
peratures and densities and we propose an interpolation 
between the different regimes with established character- 
istics, the solid and the plasma states. To take this into 
account, we have used an approach based on an inter- 
polation between current solid and plasma model follow- 
ing approach in Ref . [3(| . The relaxation rate (l/v) is 
calculated with the assumption of a screened Coulomb 
interaction among the electrons. The (e-e) relaxation 
rate rises with increasing density until it reaches a max- 
imum after which increasing screening causes the rate 
to fall. The evolution of the optical properties during 
an ultrashort laser pulse remains an open question. A 
complete description of the effect of the (e-e) collisions 
would require the complete knowledge of both the metal 
band structure and dense-plasma atomic structure, which 
makes them analytically complex and unsuitable for the 
modelling purposes considered here. In reality, it is not 
quite obvious how electron-electron collisions can induce 
a conductivity effect in the nonequilibrium solid. For a 
free electron gas, the parabolic energy-momentum rela- 
tion for which velocity and momentum are proportional 
results in neglegting (e-e) interactions contribution to ab- 
sorption. The classical approach implies that, since total 
current and total momemtum are proportional, the cur- 
rent is conserved during this type of collisions which pre- 
serves momentum. As a consequence, the conductivity 
evolution is supposed to be independent of (e-e) interac- 
tions. However, in periodic systems, as the solid sample 
in the first moments of interaction, the total momentum 
involving (e-e) interactions may be conserved via the in- 
volvement of a reciprocal lattice vector. This type of pro- 
cess will allow in consequence a variation of the electronic 
momentum and induces a change in conductivity [50l ] . A 
term v^ eq is then used to correct the conductivity dur- 
ing the nonequilibrium stage. To justify our approach, 
we expect that umklapp processes can occur at a rate 
given by a factor A with respect of the normal process 
(A < 1) [Ull. The values of this term (0.4 and 0.35 for 
aluminum and gold, respectively) have been taken from 
Refs.[53l.l54l.l55l]. being derived by low-temperature inves- 
tigations. These values will be considered constant in this 
study, although a dependence on the electronic field and 
collision rates is foreseeable. In taking into account this 
process, the slowing down of the electron flow can be ex- 
pected during (e-e) interaction and results in a change in 



the absorption coefficient when the lattice remains cold, 
during the first moments of laser irradiation. 

At near-solid density and moderate electronic temper- 
ature up to few eV, the Fermi gas undergoes collisions 
and the solid-model is given by [56[ : 



( 6_ 



/ d£(£-£ F ) 2 f(£,T e ){l-f(£,T e )} (9) 
Jo 



where f(£,T e ) is the Fermi-Dirac statistic distribution 
depending on T e and N e through the chemical potential. 



'N = -k B TJ(0,T e ) 
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(10) 
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^% where q TF = % J ^(3^N e )^ 
2ml V e e rf 



where & = m\e A /64h 3 ir 3 eQe d , m\ is the effective elec- 
tron mass 47] , A is the rate of umklapp processes, and 
Cd represents the contribution of band d electrons to the 
static dielectric function in the case of noble metals. In 
the same way as described in Ref. [HI], we have used the 
screening wave- vector qTF containing the adjustable pa- 
rameter of screening £ < 1 which was set by fitting the 
experimental results of these authors. The solid model 
is applicable up to the Fermi temperature. We assume 
that the absorption is mainly driven by the free elec- 
trons and this model allows us to calculate a straightfor- 
ward dependence in T e . This approximation provides a 
schematic view of the absorption which can be included 
in the code. For electrons at temperatures around the 
Fermi value, the phase space available for scattering in- 
creases with the electronic temperature due to the Pauli 
exclusion. Note that the increase in the (e-i) collision 
frequency due to the T e augmentation is supposed to be 
negligible with respect to the (e-e) contribution. 

For larger electronic temperatures, we use the plasma 
collision frequency given by the Spitzer formula: 



£f 
h 



£f 



-3/2 



(11) 



The approximations inherent to the plasma model are 
more acceptable at low density compared to using the 
solid model at high temperature mi- To reduce the in- 
terpolation range in the unknown region, we suppose that 
the Spitzer assumptions become suitable for T e > 2T F . 
These assumptions may be criticized because of the high 
density close to the solid one, but the choice of this di- 
rect temperature dependence is largely used to represent 
the global behavior. The nonequilibrium frequency over 
a large range of temperature is then defined by: 



vne q = Kl q (N e , T e ) for T e < T F 



'neq 



neq 



(N e ,T e ) for 7> <T e <2T F 
AN e ,T e ) forT e >2T F 



(12) 
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where v^ q = aT e 3 + &T 2 + cT e + d is a cubic interpola- 
tion which offers a realistic continuity between the two 
models. The interpolated data points are created so that 
they have a slope equal to the slope of the start or end 
segments. Note that each term (a,b,c,d) is calculated and 
tabulated as a function of N e . 




Electronic temperature (eV) 

FIG. 1: Nonequilibrium collision frequency interpolated in a 
wide range of electron temperatures in Au. Solid and plasma 
models (solid lines) are connected by a polynomial interpo- 
lation (dashed line) which allow extension in the unknown 
regime. 



As a conclusion, direct (e-e) contribution to the opti- 
cal absorption are taken into account via umklapp pro- 
cesses which are supposed to occur in the solid phase via 
a reciprocal lattice vector. This allows, in turn, a vari- 
ation of the electronic momentum, and, subsequently, 
a contribution to the optical conductivity. Concerning 
(e-i) contribution, the indirect screening of the ion po- 
tential by electronic influences is accounted for by the 
popular Spitzer dependence of the nonequilibrium colli- 
sion frequency. The equilibrium conductivity values are 
derived from theoretical considerations which match ex- 
perimental conductivity results. We suppose that the 
electric field influence is taken into account indirectly 
by the effect of the electronic energy on the collision 
rate. In this way, the nonequilibrium (T e <C T<) affects 
the Joule heating in the condensed phase in considering 
the contribution of (e-e) collisions. Although interpola- 
tion between solid and plasma regimes is questionnablc, 
we suppose that (e-e) effects prevail during the laser- 
solid interaction, yielding a dramatically modified scaling 
v neq oc T 2 (58|. A strong influence of this law is expected 
on the laser absorption, which will be discussed in the 
next sections. 



E. Interband contribution 

The absorption model was implemented for two metal- 
lic materials, gold and aluminum. For gold, no specific 
model of nonequilibrium interband contribution was im- 
plemented. We assume that the equilibrium contribution 
of the collision frequency already contains interband as- 
pects. Moreover, it is supposed that the contribution of d 
bands in the case of noble metals was taken into account 
in equation (fTUjl . 

For aluminum, the dielectric constant of the metal 
takes into account two mechanisms of laser energy de- 
position: intraband absorption in the sp band and inter- 
band absorption (parallel band). Intraband absorption is 
described through the Drude model, while the Ashcroft- 
Sturm model is used for interband contribution 59] . The 
influence of interband transition in aluminum conductiv- 
ity is then considered as cr = a D + a IB . Fisher et al. have 
reported theoretical and experimental investigations on 
aluminum absorption for femtosecond laser irradiation 
in showing the influence of the interband and intraband 
(Drude) contributions at 400 and 800 nm [36| . These re- 
sults were verified in conducting a comparable model of 
interband effects on the conductivity and similar results 
were obtained by our numerical code. 



III. EXPERIMENTAL VALIDATION OF 
THEORETICAL MODEL 

In order to validate our simulations with experimental 
results, this section is dedicated to a comparison between 
selected literature data and numerical calculations based 
on the present model of ultrafast laser absorption. The 
experimental part focuses on two demonstrations: two 
kinds of experiments were reproduced to investigate the 
optical properties of high density plasmas. Studies on 
absorption dependence in angle of incidence and inten- 
sity were performed. Moderate and high intensities were 
considered for UV ultrashort laser pulses. The results 
provide optical properties in both solid and high density 
plasma range during the laser pulse. 

Fedosejevs et al. [38| have measured the averaged re- 
flectivity corresponding of a KrF laser pulse of 250 fs at 
248 nm focused on aluminum samples as a function of the 
polarization and the angle of incidence for two intensities: 
10 14 W/cm 2 and 2.5 x 10 15 W/cm 2 . At the higher inten- 
sity, the pre-pulse is supposed to be sufficiently intense 
to produce a pre-plasma at the metal surface. In fact, 
the absorption of the main pulse is strongly modified by 
the pre-plasma occurence. In this condition we restrain 
our study to the lower intensity experiment, which corre- 
sponds to a laser fluence of 25 J/cm 2 . The experimental 
measurements were performed to demonstrate the char- 
acteristic angular dependence of electromagnetic absorp- 
tion expected for the extremely steep density gradients 
of plasmas. Our calculations can be compared with the 
experimental data since we solve the wave equation for 
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non-stationary density and temperature profiles. If the 
absorptivity decreases with the angle of incidence for the 
s-polarization, a maximum absorptivity around 62% for 
the p-polarization was observed for an angle of (54 ±3)°. 
This corresponds to the maximum energy deposited in 
the material. In the experiment, the incident beam has 
a non-desired mixed polarization and the intensity corre- 
sponding to a p-polarization is in fact composed of 93 % 
p-polarization and 7 % s-polarization. The inverse ra- 
tio corresponds to s-polarization. This polarization state 
was entered in our simulations and the calculated reflec- 
tivity is given as a function of the parameter of polariza- 
tion x p which has been fixed to be equal at 0.86 in the 
p-polarization case and —0.86 in the s-polarization case: 



100 



n a = \ 



(1 + x p )TZl + (i - Xp )n s a 



(13) 



Fig. ||2J| shows the reflectivity as a function of the an- 
gle of incidence for both polarizations. The simulated 
values were calculated with the expressions (J5)) and (fl~3)) 
with tl — 750 fs corresponding to the total pulse du- 
ration. The simulations were performed with the ex- 
pression ijSJ) for dominantly p-polarization (solid curve) 
and s-polarization (dashed curve). To show the influ- 
ence of the nonequilibrium absorption term and espe- 
cially the (e-e) collision effects, other simulation results 
are represented by setting the nonequilibrium contribu- 
tion to v equal to zero. The expression ((5J is reduced 
to the first equilibrium term which only depends on N e 
and Tj (y neq — 0). These simulations have been done 
for p-polarization (short-dotted curve) and s-polarization 
(short-dashed curve) . The reflectivity values in this later 
case are higher than the experimental ones and the sim- 
ulations that include v neq are clearly more relevant. On 
the subpicosecond timescale, the ionic temperature does 
not have enough time to increase because the energy 
transfer rate between electrons and ions has a character- 
istic time longer than tl- For the case v n eq{N e , T e ) = 0, 
the lattice remains cold during the laser pulse and the 
conductivity, depending on the electron-ion collision fre- 
quency, is close to the cold solid one. The p-polarization 
shows a minimum reflectivity at closed to 59°. This mini- 
mum which appears only for the p-polarization case, can 
be explained by the fact that for this polarization, the 
electromagnetic wave has a component colinear to the 
plasma density gradient. 

The minimal reflectivity occurs for a calculated an- 
gle of approximately 59° which corresponds roughly to 
the experimental one. Eidmann et al. have already 
compared these experimental results with their hydro- 
dynamic simulations which solve Helmholtz equation in 
a similar framework [2(j. Nevertheless, their model does 
not consider contributions of (e-e) collision in the solid 
case but the (e-ph) contribution was calculated with a 
different specific plasma model [2(| ■ In this way, the col- 
lision frequency increased quickly during the laser pulse 
and they have found an angle corresponding to the min- 
imal absorption around 70° in their calculations. They 
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FIG. 2: Simulated reflectivity as a function of angle of inci- 
dence for a 250 fs, 248 nm laser pulse for s-polarised (dashed 
line) and p-polarised (solid line) irradiation of an aluminum 
target at 10 14 W/cm 2 intensity. Minimal reflectivity is ob- 
tained for an angle of (59 ± 3)°. Simulations performed with 
Vneq = are superposed for s (short-dashed line) and p (dot- 
ted line) polarization. Experimental data points (open (s) 
and closed (p) circles with error bars) are taken from Fedose- 
jevs et al. [38j|. 



have imputed this discrepancy to the fact that they had 
not considered the depolarization in the calculation of the 
electromagnetic field. We have observed the same situa- 
tion, where the angle of maximal absorption was overesti- 
mated with a pure p-polarized pulse. In our calculation, 
this overestimation was corrected, by including a weak 
proportion of s-polarized wave in the total electromag- 
netic field. These simulations reveal a strong difference 
in the nature of the absorption of the s and p components 
of the heated matter absorption. When a density gradi- 
ent is created at the metal surface, simplest models used 
to calculate the electromagnetic fields at two interfaces or 
in the skin depth are no longer valid and it becomes nec- 
essary to solve the Helmholtz equation to calculate the 
deposited energy. Fedosejevs et al. have already indi- 
cated the importance to take into account a correct den- 
sity gradient to reproduce the experimental data. One 
dimensionnal modeling appears to be sufficient to take 
into account the effects due to incident angle and polar- 
ization. Concerning the magnitude of the reflectivity, we 
have pointed out the crucial role of the correct conduc- 
tivity by way of a total collision frequency including the 
(e-e) collision rate. Moreover, the fact that the angle of 
minimal absorption was similar to the experiment, show 
that hydrodynamic of the system is well reproduced by 
our simulations on the subpicosecond timescale. 

In order to test our model on a second type of study, we 
have investigated the laser intensity dependence of these 
two polarisations while keeping constant the angle of inci- 
dence. The dependence of the reflectivity as a function of 
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the intensity was experimentally studied by Milchberg et 
al. to determine the evolution of the aluminum resistiv- 
ity during the solichplasma transition when the change 
in density is weak [391 ]. These measurements were inter- 
preted to provide an electronic temperature dependence 
for the resistivity. This dependence was studied for an 
aluminum target with a density close to the solid one for 
an intensity evolving for four orders of magnitude, from 
5 x 10 11 to 10 14 W/cm 2 . The reflectivity was determined 
by measuring the energy after reflexion of the laser pulse 
which hits the target at an incident angle of 45°. The 
target were prepared by depositing a 400 A aluminum 
film on a glass substrates. The authors stated that for 
thicker film, the reflectivity remains constant. 




Intensity (W/cm 2 ) 

FIG. 3: Simulated reflectivity as a function of intensity in an 
aluminum sample for s (dashed line) and p (solid line) polar- 
ization. Simulations performed with v neq — are superposed 
for s (short-dashed line) and p (dotted line) polarization. Ex- 
perimental data points (black dots) are taken from Milchberg 
et al. [39] . 



Fig. ([3]) shows the experimental reflectivity associated 
to the 400 fs (FWHM) s- and p-polarized incident pulses 
for fluences from 20 mJ/cm 2 to 2 kJ/cm 2 at X a = 308 
nm. A decrease of reflectivity with intensity is observed. 
The challenge is to relate this decrease to an excitation 
model. The figure also shows the calculated result in two 
situations [v neq = 0, v neq ^ 0). Our simulations are in 
good agreement with these experimental measurements 
when the (e-e) collisions are considered. Numerical cal- 
culations provide reflectivity evolutions based on the ex- 
pression ([S]) where the upper limit of the integration time 
corresponds to the total pulse duration tl = 1.2 ps. A 
drop of reflectivity for an intensity range of 10 12 — 10 15 
W/cm 2 is observed. In the case v neq = 0, the dependen- 
cies in Tj and iV e on the (e-i) collisional processes are 
not sufficient to reproduce the experimental data and 
the agreement is poor. The drop in reflectivity occurs 
for higher intensities for both polarizations when the in- 
crease in ionic temperature starts to be important and 



to drive the absorption mechanism. This again shows 
the role of considered (e-e) collision in influencing the 
transient optical properties. 

Three regimes are clearly distinguishable when observ- 
ing reflectivity as a function of the laser intensity. The 
reflectivity can be written as a function of the refractive 
index n as TZ = |(1 — n)/(l + n)\ 2 , where, according to 
the Drude expression: 



n {w a jv) = 1 + 



iui a jv 



(14) 



where 



(N e e /m' e eo) is the plasma frequency. The 



above equation used in the case of normal incidence of 
laser radiation, and the reflectivity at other incident an- 
gles for p- and s-polarized radiation can be calculated 
based on the Fresnel formulas. These considerations do 
not bring new informations for a quantitative discussion 
but allows a qualitative one. If 1 < uj a /v < u> p /v with 
uj a = 6.12 x 10 15 s _1 , then h 2 ~ 1 + iL0 2 /uj a v and lZ a 
decreases with is, corresponding to the first regimes. The 
first one, indicated by the label (a) is observable at inten- 
sities lower than few 10 12 W/cm 2 . This regime, caracter- 
ized by a reflectivity which decreases slowly, correspond- 
ing to an electronic temperature which is not sufficiently 
high to induce high collision frequency (y <§; u> a ). More- 
over, the energy exchange between electrons and ions is 
low and the matter remains in a solid state during the 
pulse duration. The TZ a decrease is related to the rise of 
(e-i) collision frequency, which depends almost linearly 
on the ionic temperature. Moreover, the decrease in den- 
sity is weak in this regime and the reflectivity does not 
depend on the electronic density through v and ui p . The 
second regime, marked by the label (b) on the Fig.®, is 
reached for intensities between 10 12 and 10 14 W/cm 2 . In 
this case, v ~ tv a and the drop of lZ a can be related to 
a strong T e increase which leads to a frequency collision 
augmentation. In the assumption of the solid-model, lZ a 
decreases as v neq cx T 2 , which corresponds to the regime 
put forward by Milchberg et al. in Rcf . [39] . Note that the 
decrease in density starts to play a significant role in this 
regime as well and can induce a reduction of the reflec- 
tivity. Finally, for the highest intensities, the reflectivity 
remains constant or increases. This regime is more diffi- 
cult to reproduce with simulations because, beyond 10 14 
W/cm 2 the experimental pre-pulse is sufficient to pro- 
duce vapor or plasma before the absorption of the main 
pulse. In this case marked by the label (c), the electronic 
critical density reached in front of the aluminum sample 
yields a strong modification of the absorption proper- 
ties. For the p-polarization case, 1 <C lo p /v < Lo a /v and 
lZ a ~ (ujp/2uj a ) 2 becomes low, depending on the elec- 
tronic density. 

The good agreement between the experimental and 
theoretical reflection factor shows the possibilities of us- 
ing the presented model to reproduce and analyse the 
optical properties during the plasma formation induced 
by a laser pulse. This was achieved by properly estimat- 
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ing the collision frequency and other optical and ther- 
modynamical parameters in noncquilibrium conditions. 
These results indicate that the numerical code is able to 
simulate the absorption of the incident electromagnetic 
wave in electron-ion noncquilibrium regimes. This pro- 
vides opportunities to caracterize the response of a mate- 
rial for timescales lower than the electron-ion relaxation 
time. We have studied reflectivities induced on a thick 
metal sample, where the heated matter at the surface 
was composed by a temporal sequence of several layers 
of different thermodynamic states. The next section is 
dedicated to specific experiments which allow to follow 
a single thermodynamic state on ultrashort timescales 
ensuring as well temporal resolution. Our simulations 
are suitable to provide precious informations on the col- 
lisional processes and state of the system during its re- 
laxation. 



IV. APPLICATION TO PUMP-PROBE 
DIAGNOSTIC 

Optical reflectivity and transmission measurements 
provide valuable insights into studying ultrafast phe- 
nomena in condensed matter and plasmas physics [1 SI ] - 
In particular, they enable observations of phase transi- 
tions and provide information about the transport prop- 
erties of dense plasmas. If the resolution is short enough, 
pump-probe experiments are relevant to characterize the 
properties of the material during relaxation processes. 
The use of ultrashort laser pulse provide a sub-picosecond 
resolution of the evolution after an ultrafast laser event, 
the pump excitation. The subsequent dynamics is then 
governed by the interactions among the elementary exci- 
tations and is monitored by using a second optical pulse, 
the probe. Reflectivity measurements strongly depend 
on plasma density and collision frequency. In these con- 
ditions, numerical investigations can provide the corre- 
lation between matter properties and relaxation process. 
In this section, we compare experimental results and sim- 
ulations performed in similar conditions which permit us 
to split up the sequence of electron-electron, electron- 
phonon and hydrodynamics relaxations. 



A. Strong pulse and unperturbing probe 

Our assumptions are firstly tested against the data 
points corresponding to pump-probe reflectivity and 
transmission experiment performed by Widmann et al. 
in Ref. §. In this case, the target is a gold foil with a 
thickness of about 28 nm which has been irradiated with 
a A Q = 400 nm, 150 fs (FWHM) laser pulse, at normal 
incidence. After focusing, the intensity of the pump pulse 
is about 10 13 W/cm 2 on the metal surface. The reflec- 
tivity and the transmission of the initial solid and the 
produced plasma are analyzed with a s-polarized probe 
pulse at \p = 800 nm with a similar duration, an in- 



cident angle of 45°, and a time delay varied from to 
15 ps after the pump pulse. In the experiment, reflected 
and transmitted energies were measured with photodi- 
odes and compared to the incident energy to deduce an 
instantaneous reflectivity and transmission of the mat- 
ter. In a similar way as in the previous section, simu- 
lated reflectivity TZg(St) and transmission ^ s p(5t) were 
calculated with the expressions ([5]) and ([6]) where r was 
taken to be equal to a variable timestep 5t ^ 0.1 fs. 
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FIG. 4: Numerical evolution (curves) of the reflectivity and 
the transmission of a 800 nm pulse probing a 28 nm gold foil 
thickness, for different simulation configurations. The solid 
line (a) corresponds to the best-fitted complete simulation. 
The dash-dotted line (b) is a simulation without nonequilib- 
rium contribution in the absorption. The dashed line (c) and 
the short-dotted line (d) correspond to simulations performed 
with a lower (7 = 7o/4) and a higher (7 = 70 x 2) e-ph cou- 
pling term, respectively. The hydrodynamics processes have 
been switched off in the simulation shown by the short-dashed 
line (e). Experimental data points (black diamonds) are taken 
from Widmann et al 

At the beginning of the irradiation and at the probe 
frequency, solid gold shows a reflectivity of 93%. Follow- 
ing the irradiation, three regimes are clearly distinct on 
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the Fig. ([4]). The first one, shortly after the initial mo- 
ment, is characterized by a swift decrease in both reflec- 
tivity and transmission during the first picosecond. The 
second one is marked by a plateau between 1-7 ps and 
corresponds to a regime where optical properties seems to 
be frozen because reflectivity and transmission are almost 
constant. Finally, a change appears after 7 picoseconds 
and a last regime starts, characterized by a decrease in 
reflectivity and an increase in transmission. To explain 
these behaviors, we have performed five simulations with 
different configurations to identify the effect of the sup- 
posed main processes. 

The first simulation labelled by (a) on the Fig. ((4J is a 
complete simulation including all the processes described 
in the section II and is used as the reference case to inter- 
pret the other calculations. The dash-dotted line (b) is 
a simulation without noncquilibrium contribution in the 
collision frequency. In this contribution has 

been skipped from the equation J5]) and no (e-e) collision 
mechanism is taken into account in the bremsstrahlung 
absorption calculation. The expected difference has to 
take place during the nonequilibrium process, before 
electron-phonon thermalization. The calculated reflec- 
tivity is strongly different in the case (b) since the con- 
ductivity parameter only depends on the density and on 
the ionic temperature which evolves slowly compared to 
the electronic temperature variation. In contrast, curve 
(a) shows a drop in the reflectivity and reaches a mini- 
mum value at the end of the laser pump pulse, when the 
electronic temperature is maximal. At this time, SJie^] 
and the subsequent quantity of absorbed energy are max- 
imal. The transmission curve displays a less strong dif- 
ference for these two cases because IZp diminution and 
absorption increase occur simultaneously. For the case 
(a), as electrons transfer their energy to the ions, v de- 
creases and Tip rises and reaches a maximal value at 4 
ps. From this study, we can infer that the observed ab- 
sorption in ultrafast timescales is mainly due to (e-e) col- 
lision. Thus, during the electromagnetic excitation, the 
Coulomb interaction between electrons drives the Joule 
heating process. 

The curves (c) and (d) include nonequilibrium absorp- 
tion mechanism but differ from (a) in ( e-ph) energy trans- 
fer rate. In fact, the case (c) corresponds to a case where 
the usual (e-ph) coupling term 7 = 4 x 10 16 W K _1 m~ 3 
has been divided by a factor 4 to slow down the ionic 
temperature increase and to increase the nonequilibrium 
lifetime. Compared to the experimental data and the (a) 
reference case, the break-up with the third stage at 7 ps 
is smoothened and the region where the reflectivity re- 
mains constant is more pronounced. In contrast, for the 
(d) case, the (e-ph) coupling term 70 was multiplied five 
times to accelerate the ionic temperature increase and to 
decrease the nonequilibrium lifetime. In this case, the 
length of the plateau (which is visible in both reflectivity 
and transmission) decreases and the equilibration time 
seems to be reduced. We deduce that the plateau can be 
attributed to the temperature equilibration mechanism 



between electrons and ions and the duration of this stage 
corresponds to the (e-ph) relaxation time. After this, the 
hydrodynamics commences. 

Finally, the fifth round (e) was performed by keeping 
7 = 70 as in (a) but the hydrodynamics processes have 
been stopped. To do this, the pressure p and the density 
p equal the standard pressure po = 10~ 5 Pa and the 
solid density po — 19-3 g cm~ 3 , respectively, in the entire 
gold film. In this non-physical case, no pressure gradient 
can cause motion, volume increase, surface expansion, or 
shock processes in the solid. The behavior of the optical 
properties is similar for (e) and (a) during the 3-4 first 
picoseconds and, after this time, the reflectivity and the 
transmission split and follow very different evolutions. In 
the isochoric case, the reflectivity increases continuously 
as the energy is transferred from the electrons to the ions. 
In fact, the total collision frequency v decreases because 
the augmentation in (e-i) collision is not compensated 
by the (e-e) collision decrease. The transmission does 
not evolve anymore after the end of the nonequilibrium 
process since the change in reflectivity is compensated by 
a weaker absorption of the probe. In the hydrodynamic 
case (a), the reflectivity decreases because of the density 
decrease which is related to the plasma formation at the 
surface. In fact, according to the Spitzer formula and 
the Drude model for v < cup, IZp decreases with the 
ionic temperature increase due to the increase in (e-i) 
collision frequency and IZp decreases when the electron 
density drops. The change in 073 yields a greater skin 
depth for the probe, leading to increases in optical path 
length. As the plasma expands, the temperature and the 
density falls and the free electrons recombine. 

The experimental measurements were performed on 
30 ± 2 nm gold sample. The magnitude of absorption 
strongly depends on the material thickness. To estimate 
the influence on our simulations, Fig. [5] shows a depen- 
dence of the transient optical properties in the material 
thickness. In our simulation, an increase of energy of 
about 4.8 x 10 6 J/Kg notified by Ref.[|[ corresponds to 
an input fluence of 0.5 Jem -2 for the pump pulse inci- 
dent on a sample of 28 nm thickness. The behavior of 
reflectivity and transmission were investigated for three 
thicknesses: 28, 30 and 32 nm at the same fluence. These 
three simulations show that the general evolution of the 
optical properties are identical but differs in magnitude. 
The reflectivity increases whereas the transmission de- 
creases with the film thickness. This simple result can 
give an explanation to the discrepancy between the ex- 
perimental and simulated values in providing a kind of 
error bar to the numerical experiment data. 

This numerical study emphasizes the role of the (e-e) 
umklapp collision mechanism on the absorption of the 
laser pulse. The (e-e) interaction influence the mate- 
rial optical reponse during a short time after the pump 
pulse. The following (e-i) equilibration and the hydro- 
dynamic stages have been identified as well. Due to the 
number of approximations in the optical model, the sim- 
ulation does not perfectly match the experimental data. 
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FIG. 5: Time evolution of the reflectivity and the transmis- 
sion of a 800 nm probe pulse for different thicknesses of a 
gold sample pre-excited by a 400 nm laser pulse. Both be- 
haviors exhibit energy relaxation stage and solid to plasma 
transformation induced by a 400 nm, 0.5 Jem -2 pump pulse. 
Experimental data points (black diamonds) are taken from 
Widmann et al [5J. Generation of nonequilibrium, relaxation 
to equilibrium and solid to plasma transition are visible. 



B. Experiments under near-threshold laser 
excitation on solid samples 

The final discussion relates to a class of pump-probe 
experiments which focus on particle emssion. We pro- 
pose a direct application of a double-pulse sequence in 
studying the consequence of transient optical properties 
on energy coupling efficiency as a function of delay. Ex- 
periments based on ion signal emission as a function of 
the time delay were performed on several metals [60j. 
The experiments map in time the energy deposition as 
a function of the delay between identical sub-threshold 
pulses, the total energy being above the threshold. Here, 
the threshold is defined as ion emission threshold. This 
way a quantitative picture of the transient behavior of 
the energy coupled into the sample emerges. It is sup- 
posed that the energy density stored in the metal surface 
is strongly correlated on the ion emission (25|. A fast 
decaying ion signal that drops to zero in the first hun- 
dreds of femtosecond [6(| was followed by a large peak of 
emission in the picosecond range [25|, [60, [6l[ . 
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Nevertheless, the different phases of the matter evolution 
and the subsequent changes in the behavior of the opti- 
cal properties are quite well reproduced. This numerical 
study devoted to reflect the experimental conditions in- 
volving laser-metal irradiation for both pump and probe 
pulses is an efficient tool to discriminate collisionnal pro- 
cesses. A direct interest of a low energy probe pulse is 
the possibility of carrying out tests of optical absorption 
of the material. This study could be extended to perform 
temporal optimisation of laser absorption in determining 
the time when the matter state is the more prepared and 
able to respond to the electromagnetic excitation. 



FIG. 6: Calculated specific absorbed energy in the skin depth 
of an aluminum target irradiated by a double pulse (2 x 0.4 
Jem -2 ) and a gold target irradiated by a double pulse (2 x 1.9 
Jcm~ 2 ). Simulations were performed as a function of time 
delay. 



To discuss the experimental results related to the time 
varying ion signal, we have calculated the specific ab- 
sorbed energy inside the skin depth in a metallic sam- 
ple. We have chosen Al and Au as examples. The skin 
depth is defined as a time-dependent quantity, depend- 
ing on the imaginary index which evolves with space and 
time. The laser irradiation is composed by two numer- 
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ical identical pulses with a varying delay. We make the 
assumption that, to a certain extent, the energy density 
in the superficial layers determines the efficiency of phase 
transformation into a gas state, and, therefore, particle 
emission [25[. The total specific absorbed energy (E sa ) 
of a 800 nm, 150 fs double-pulse is then calculated as 
a function of time separation between the two pulses. 
We discuss two cases, Al and Au for two reasons. The 
first reason is related to the different values of (e-ph) 
coupling, Al showing an almost 10 times higher electron- 
phonon interaction strength. Secondly, Al shows a par- 
ticular optical characteristic, with no significant change 
between solid and liquid phases at 800 nm [62j. The nu- 
merical experiment are carried away with two identical 
sub-threshold pulses of 150 fs and 800 nm wavelength. 
The threshold is defined numerically by the onset of the 
gas-phase at 0.45 Jcm~ 2 for Al and 2 Jem" 2 for Au. The 
results are presented in Fig. ©. For both material, the 
following particularity were noted. 

Several characteristic stages are noticed with various 
peaks in the absorbed energy density. The first one ap- 
pears during the first picosecond and shows a decay in 
the amount of energy stored E sa , due to the decrease of 
3?e[^]. This decay is faster for Al than for Au. Note that 
an important result not visible here is that simulations 
with a single pulse and a nonperturbing probe have re- 
vealed that the absorption increases strongly during the 
first pulse. This increase can not be observed with an 
identical double-pulse experiment due to a convoluting 
effect. This increase is a consequence of the instanta- 
neous T e increase which yields a strong augmentation 
of the collision frequency. As T e reaches its maximum, 
the absorption is maximal and the observed minimum is 
almost synchronized with the corresponding time. The 
electron-phonon energy transfer is responsible for the T e 
decrease. As a consequence, T e can be maximal at the 
end of the first pulse if the electron-phonon coupling con- 
stant is very low. In the cases presented here, the absorp- 
tion is maximal at about the half of the pulse duration, 
the corresponding time delay is close to the axis origin 
and a fast decay occurs for short delay times. This vari- 
ation is related to the T e decrease which induces a drop 
of the collision frequency v and a diminution of 5R e [erg] . 
According to eq. in the range up/v > 1, Ke[crg] an d 
the associated Joule heating process decreases as v de- 
creases. Moreover, a simulation performed, with setting 
v n eq = 0, has shown that there is no decay for short time 
delay in this case. The electron-phonon energy transfer 
rate is correlated to the slope of the E sa curve at short 
times. This slope is steeper for Al than for Au and the 
decrease ends at 1.5 ps for Al and 6 ps for Au. A liquid 
phase is reached around 2 ps for Al and 7 ps for Au. The 
effect of the solid-liquid transition on the optical proper- 
ties, which is expected to be more important for Au than 
for Al, is negligible because the nonequilibrium mecha- 
nism of absorption, especially (e-e) umklapp collisions, 
dominates. During this time, the ionic temperature Ti 
increases and the v el increase results in an augmentation 



of SRe^], balancing the effect of T e decrease. At this 
stage, the optical properties are dominated by the ionic 
heating and Ti drives the absorption mechanism. After a 
single pulse, the thermalisation is reached around 15 ps 
for Al and 50 ps for Au. Ti begins to decrease while p 
augments, prolonging the absorption increase. Note that 
solidification process begins to occur at the rear side of 
the liquid layer. At 35 ps for Al and 75 ps for Au, solid- 
ification reaches the free surface of the liquid layer and 
p remains stable at solid density. The temperature de- 
crease induces a drop of E sa for longer time delays. 

These simulations have shown that a fast decay in the 
absorbed energy can be imputed to the electronic tem- 
perature decrease and associated (e-e) effects during the 
nonequilibrium stage. Due to a smoother lattice heating, 
the fast decay is slower for Au than for Al. A recov- 
ery of the energy coupling efficiency on a tens picosecond 
scale is also observed for both materials. We have shown 
that the fact that the enhancement of optical absorption 
on picosecond scale is due to an elevation of the ionic 
temperature. This enhancement happens earlier for Al 
because of a high electron-phonon coupling constant. Re- 
condensation of matter into a high density phase occurs 
after several tens of ps and result in a E sa drop. A similar 
situation was observed in experiments. As a conclusion, 
umklapp electron-electron collisions in the solid state de- 
termine the optical absorption behavior when the fluence 
is close to the threshold value. 



V. CONCLUSION 

Collisional absorption by inverse bremsstrahlung plays 
a major role in laser-matter interactions. The full de- 
scription of inverse bremsstrahlung absorption requires 
knowledge of the electron-photon absorption mecha- 
nism and the associated umklapp electron-electron and 
electron-ion collision rates. We have conducted a sim- 
plified modelling of the optical properties to perform re- 
alistic simulations of several experiments present in the 
literature. We have used a relaxation-time model to ap- 
proximate the two-body scattering rate. An important 
parameter determining the modification of laser-light ab- 
sorption is the conductivity parameter and we have pro- 
posed a description of its evolution during the material 
excitation. We used published results to test the hy- 
pothesis of an electronic contribution in an unified frame, 
covering various excitation stages, against its ability to 
explain experimental observations. In bringing together 
several accepted models in solid and plasma range, a con- 
ductivity dependence of macroscopic properties such as 
electronic temperature (T e ), ionic temperature (Ti) and 
electronic density (N e ) was employed. This study is in- 
dicative to the potential role of electron-electron umklapp 
processes in optical absorption. An accurate evaluation 
of the umklapp rate would require complex kinetic ap- 
proaches. 

The numerical diagnostic that we have developped pro- 



13 



vide the correlation between the sequence of states of 
matter and the optical response. The information, re- 
solved in time, give insights into the relaxation mech- 
anisms which drive the phase transitions and open new 
opportunities to monitor the excitation of a metal in con- 
troling both its optical and thermal evolution. Several 
regimes of irradiation, including above and subthreshold 
regimes were investigated. Clear modification of elec- 
tronic and hydrodynamic effects was observed. Beyond 
the retrieval of the optical properties and the discussion 



of associated thermodynamic states, a qualitative idea on 
the transient energy coupling properties was given. 
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